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Abstract 

New Hamiltonian formalism based on the theory of conjugate curvilinear coor- 
dinate nets is established. All formulas are "mirrored" to corresponding formulas 
in Hamiltonian formalism constructed by BA. Dubrovin and S.P. Novikov (in a 
flat case) and E.V. Ferapontov (in a non-flat case). In the "mirrored-flat" case 
Lagrangian formulation is found. Multi-Hamiltonian examples are presented. In 
particular Egorov's case, generalizations of local Nutku-Olver's Hamiltonian struc- 
ture and corresponding Sheftel-Teshukov's recursion operator are presented. An 
infinite number of local Hamiltonian structures of all odd orders is found. 
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1 Introduction 



The modern theory of integrable hydrodynamic type systems 
of first-order PDE's (see e.g. [50], [51]) 



i.e. quasilinear systems 



v\{u) 



V{[U) 



v-[u) 



( ul \ 



V < / 



(1) 



where u l = u l (x,t), i = 1, • • • ,n are the field variables (unknown functions), x and t are 
one-dimensional space and time variables — started with the paper [14] where a general 
Hamiltonian formalism for such systems was proposed. Physical examples of (1) with a 
"natural" Hamiltonian structure are numerous (see e.g. [3], [15], [30], [31]) for physical 
examples of systems with Hamiltonian structures originating from infinite-dimensional Lie 
algebras naturally related to the respective systems; such Hamiltonian structures are a 
particular case of the more general Hamiltonian structure described by B.A. Dubrovin and 
S.P. Novikov. Their generalization allows us to cast into the Hamiltonian form another 
vast class of hydrodynamic type systems (1), derived from integrable nonlinear equations 
via the so called "averaging" procedure, which gives a system (usually of type (1)) de- 
scribing slowly varying ("modulated") quasiperiodic solutions (see e.g. [9], [14], [18], 
[24], [33], [35] [56], [57]). These Whitham equations possess the Hamiltonian structure 
of the Dubrovin-Novikov type. Some particular cases of this Hamiltonian structure for 
averaged equations were studied in [29]. 

This method may be applied to a wide range of "averaged" integrable equations which 
possess alongside with the Hamiltonian structure of the Dubrovin-Novikov type the re- 
markable property of diagonalizability: after a suitable change of the field variables r l (u) 
1) one obtains the hydrodynamic type system 



in 



v r r„ 



l,2,..,iV. 



(2) 



(no summation over repeated indices hereafter!) with diagonal matrix v l (r)5j. These new 
variables r % are called Riemann invariants (in the Whitham theory they are branch points 
of corresponding Riemann surfaces, see e.g. [24], [48]). For n = 2 every hyperbolic 
system (1) may be diagonalized; for n > 2 this is not true in general, a special criterion 
[28] shall be applied to check diagonalizability of (1). 
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On the other hand the methods of [56] (see also [57]) lead to an unexpected link 
between the theory of Hamiltonian diagonalizable hydrodynamic type systems and a clas- 
sical object of local differential geometry intensively studied in the end of 19th and be- 
ginning of 20th century — orthogonal curvilinear coordinates in the flat Euclidean space 
R n . Practically all "physical" entities (conservation laws, symmetries etc.) have their 
counterparts in the theory of orthogonal curvilinear coordinates; amazingly enough the 
basic formulas relating conservation laws, symmetries and the Hamiltonian structure of 
integrable systems of hydrodynamic type may be found in [10]! Further investigation 
discovered a deep relation between other types of integrable nonlinear systems studied 
in the modern theory of integrable physical equations and classical problems studied in 
[5, 6, 10, 11, 12, 17, 27] and other papers at the beginning of previous century 

In this paper we expose this remarkable link as well as some recently discovered deeper 
relations between the classical differential geometry flourished at the end of the 19th cen- 
tury and the modern theory of integrable nonlinear PDEs appearing in different applica- 
tions to mathematical physics. 

The theory of semi-Hamiltonian hydrodynamic type systems is a most developed part 
of integrable PDE systems (see [14], [56], [57]). A one of such reasons is that this theory 
is based on classical differential geometry, and more precisely, on the theory of conjugate 
curvilinear coordinate nets (see [10]). For instance, the local Hamiltonian formalism 
associated with the differential-geometric Poisson brackets of the first order is connected 
with theory of orthogonal curvilinear coordinate nets (see [13], [32], [56], [57], [59]). 

This paper is devoted to an alternative construction based on so-called anti-flatness 
property. An existence and, moreover, an explicit construction of Hamiltonian structures 
play important role in integrability of hydrodynamic type systems, because a general 
solution of semi-Hamiltonian hydrodynamic type system depends on general solution of a 
linear problems describing conservation laws and commuting flows, which are related by 
the Hamiltonian structure. 

In general case the problem is following: it is necessary to describe all possible Hamilto- 
nian formalisms associated with integrable hydrodynamic type systems or with conjugate 
curvilinear coordinate nets. The local Hamiltonian formalism determined by a differential 
operator of the first order is constructed in ([14] (see also [56], [57]). The nonlocal Hamil- 
tonian formalism of the first order associated with a flat normal bundle is constructed 
in [19]. This paper is devoted to the nonlocal Hamiltonian formalism associated with 
the local Lagrangian representation of corresponding hydrodynamic type systems. Sup- 
pose a given hydrodynamic type system possesses the local Hamiltonian structure of the 
Dubrovin-Novikov type (see [14]) and the above mentioned Lagrangian representation, 
simultaneously. It means, that this hydrodynamic type system has local Hamiltonian and 
local symplectic structures. A most amazing consequence is that the corresponding hy- 
drodynamic type system has infinitely many local Hamiltonian structures of all odd orders 
(as well as infinitely many local symplectic structures and corresponding local Lagrangian 
representations). 

This paper is organized as follows. In Section 2, the relationship between local Hamil- 
tonian formalism of the Dubrovin-Novikov type and orthogonal curvilinear coordinate 
nets is briefly described. In Section 3, the first "puzzle" hidden in the paper [42] written 
by Yavuz Nutku and Peter Olver is formulated. In Section 4, integrable hydrodynamic 
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type systems determined by a special local Lagrangian representation are found. The 
existence of such Lagrangians is based on the concept "anti-flatness" , which is established 
in this paper. In Section 5, simultaneously, hydrodynamic type systems associated with 
a local Hamiltonian structure and with above mentioned local Lagrangian representation 
are considered. Local Hamiltonian structures of all odd orders are described. In Sec- 
tion 6, hydrodynamic type systems possessing anti-flat multi-Hamiltonian structures are 
discussed. In the Section 7, the above construction is restricted on the symmetric case. 
Statement: suppose the Egorov hydrodynamic type system has just one known local 
Hamiltonian structure of the Dubrovin-Novikov type (or vise versa has just one known 
local Lagrangian representation), then infinitely many other local Hamiltonian structures 
of all odd orders can be constructed. All of them can be expressed via corresponding 
solutions of the WDVV equation (see [13], [34], [38], [46]). In Section 8, the second 
"puzzle" hidden in the paper [42] written by Yavuz Nutku and Peter Olver is formulated. 
A natural interpretation via local Hamiltonian structures of the Dubrovin-Novikov type 
is given. In Section 9, a generalization of above construction on a nonlocal case associated 
with surfaces with a flat normal bundle is briefly discussed. In Section 10, the above local 
Lagrangian formulation is extended on higher order commuting flows of hydrodynamic 
type systems. 

2 Flat case 

The theory of conjugate curvilinear coordinate nets (see [10]) is based on nonlinear PDE 

system 

di(3 jk = PjiPik, i^j^k, (3) 

where rotation coefficients (5 ik are functions of N independent variables r k (which are 
called principal curvature coordinates, see [54]). A general solution of this system is 
parameterized by N(N—l)/2 functions of two variables. The above nonlinear PDE system 
is a consequence of the compatibility conditions di(d k Hj) = d k (diHj), di(d k ipj) = d k (diijjj), 
i/j/ k, where Lame coefficients Hj and adjoint Lame coefficients tpj are solutions of 
a linear PDE system and an adjoint linear PDE system, respectively 

diH k = p ik H u dktpi = p ik ip k , % ^ k, (4) 

which have general solutions parameterized by N arbitrary functions of a single variable 
for any given rotation coefficients (3 ik satisfying (3). In general case a relationship between 
these linear PDE systems is unknown. 

A link connecting conjugate curvilinear coordinate nets and integrable hydrodynamic 
type systems is following. 

• Let Hi and Hi be two particular solutions of the first linear PDE system (4). 

• Let field variables r k be functions of two independent variables x and t. 

• Let construct a hydrodynamic type system written via Riemann invariants r k 



Let field variables r k be functions of three independent variables x, t and r, simul- 
taneously. Then the above hydrodynamic type system possesses a commuting flow 
(i.e. (r*) t = (rj) T ) 

r* = §4 (6) 
parameterized by N arbitrary functions of a single variable (see (4) and [56], [57]). 
Let introduce the function h such that 

dih = ipA- (7) 

Then the above hydrodynamic type system has a conservation law (written in the 
potential form) 

d£ = hdx + gdt + fdr, 

where 

= ifrfii, dif = ip^. 

Thus, these hydrodynamic type systems possess infinitely many conservation laws 
parameterized by N arbitrary functions of a single variable. 

A general solution of the hydrodynamic type system (4) is given (in an implicit 
form) by the generalized hodograph method (see [56]) 

xH, + tHi = H. 



The first problem in the theory of hydrodynamic type systems (5) is a description of 
solutions of the nonlinear PDE system (3), the second problem is a description of solutions 
of the linear PDE systems with variable coefficients (4). These problems are very compli- 
cated in general. By this reason, conjugate curvilinear coordinate nets can be separated 
on some sub-classes due to some criterion (features, characteristics etc.), then corre- 
sponding conjugate curvilinear coordinate nets can be integrated within corresponding 
sub-classes. 

The approach developed by B.A. Dubrovin, S.P. Novikov (see [14]) and S.P. Tsarev 
(see [56]) is based on a local Hamiltonian formalism for hydrodynamic type systems (5) 
connected with the theory of orthogonal curvilinear coordinate nets (see [56], [57]), 
which is a sub-class of the theory of conjugate curvilinear coordinate nets. 

• Let establish a link of the first order between solutions of the linear and adjoint 
linear PDE systems (4) 

Hi = diipi + ^2l3 mi ip m . (8) 

• Then the flatness condition is given by (see [10]) 

diPik + d kPki + Yl PniiPmk = 0, i ± k (9) 
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which can be obtained by the substitution (8) in (4). Thus, the nonlinear PDE 
system (3) and (9) is a consequence of compatibility conditions of the linear PDE 
system (4) and (8). It means that orthogonal curvilinear coordinate nets are de- 
scribed by the nonlinear PDE system (3) and (9). 

If l.h.s. of (8) is vanished, then N particular solutions Vv 7 )j are determined by the 
linear ODE systems (for each fixed index i; see (4) and (8)) 

These linear ODE systems have N(N + l)/2 constraints ("first integrals") 

9fh = Yl ^G8)T*^(7)m = const • ( 10 ) 
This is non-degenerate matrix. Then = g^ip^. 
Let a 7 be flat coordinates, i.e. <9,a 7 = ^jy^Hi (see (7)) and 

<9a 7 Hi ' 



'IE 



Then (5) under a point transformation r % = r J (a) can be written in the conservative 
form 

where h(a) is a Hamiltonian density, the momentum density h(a) is a quadratic 
expression with respect to flat coordinates a 7 

h = \g^a?a>. (13) 

It means (see [49]), that the above hydrodynamic type system has two extra con- 
servation laws (of the momentum and of the energy, respectively) associated with 
the above local Hamiltonian structure 

Then the above conservation law densities h and h can be expressed via Riemann 
invariants r k (in quadratures) 

dih = d t \f^^) =^H t , di ^L = ^ m . Hh 



where (see (8)) 

Hi = di^i + ^Pmi^m, Hi = difa + ^/3 m ^ m - (14) 

Corollary: The momentum density (see the first above equation) is given by (cf. 
(13)) 

h = \^l- ( 15 ) 

It means, that the local Poisson bracket of the Dubrovin-Novikov type in flat coor- 
dinates is given by 

{a f> (x),a?tf)} = g (hf 6\x-x'); 

in Riemann invariants 

{r i (x),r j (x')} = A ij 6(x-x'), 

where 

A " = W d *W> = WW^kd - P ik r k x ). (16) 

Hi Hi tiiH k 

Thus, an arbitrary commuting flow (6) also can be written in the same local canon- 
ical Hamiltonian form 

where 

and solutions and Hi of linear PDE systems (4) are related by (8). 

The relationship of the first order between the linear problems (4) is given by (8). 
Taking into account another identity 

an inverse (nonlocal) relationship between the linear problems (4) is given by 



da^ 3 ' 
Indeed (see (8)), 
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3 First Nutku— Olver's "puzzle" 

The ideal gas dynamic system (see e.g. [41], [42]) 



u 2 



p t = (pu) x , u t = ^- + pf'(p) - f'(p)j (19) 

has only one local Hamiltonian structure of the Dubrovin-Novikov type 

p\ _ ( D x \(Qfr 



u I V D r 



s 



where the Hamiltonian and other lower functionals are given by 

,2 



^ r + pf>(p)-2f(p) 



dx, H3 = / pudx, H2 = / pdx, Hi = / udx 



However, Ya. Nutku and P. Olver (see [42]) found an absolutely new local Hamiltonian 
structure determined by the third order purely differential operator 

B = W n (20) 



, -Dr 

where the recursion operator R is a purely differential operator of the first order 

R = D x (W x y 1 (21) 

and the matrix W is given by 

W=( 2, P \- (22) 



J"(p) u 

This paper is devoted to an investigation of this phenomenon. Indeed, this local 
Hamiltonian structure of the third order has clear pure differential-geometric meaning. 
Moreover, a corresponding differential-geometric construction is generalized on iV com- 
ponent hydrodynamic type systems. 

So, we have the set of questions: 

• How to explain the ORIGIN of this local Hamiltonian structure? 

• How to generalize this local Hamiltonian structure on N component case? 

• How to recognize that a given hydrodynamic type system possesses such 
a local Hamiltonian structure? 

• How to construct such a local Hamiltonian structure? 

• How many such local Hamiltonian structures a given hydrodynamic type 
system can possess? 



S 



• Is it any relationship between these local Hamiltonian structures and 
local Hamiltonian structures of the Dubrovin— Novikov type? 

The key idea used for opening this puzzle is hidden in the local Lagrangian represen- 
tation for this hydrodynamic type system (see [43]) 



S 



1 p x u t -u x p t 



2ul-f"'{p)pl 

4 Local Lagrangian representations 

The Lagrangian 

S =/ b4 (r,r„r„...) f f- fc (r,r.,r„...)]^ 



dxdt. (23) 



determines the Euler-Lagrange equations 
in the Hamiltonian form (see [40]) 

r i = K ij — 

where the Hamiltonian H = J h(r,r x ,r xx , ...)dx, and the Hamiltonian operator K l i is an 
inverse operator to the symplectic operator 

4 = • (24) 

Let us restrict our further consideration on the case 

S =/to(r,r. W -Mr,r.,r.,...)]dx*. 



Then corresponding Euler-Lagrange equations can be written in the form 
where a corresponding symplectic operator is given by (see (24)) 



dr l dr k \ dr% dr 

If we choose 
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then corresponding components of the above symplectic operator are given by (cf. (16)) 



m.. = Ei d Ei 

lvl ll 1 i 



M ik \k^i — HiH k 



hi 



where 







ik 



djH k 



k^i. 



Theorem 1: The Hamiltonian operator 



n — e w^r x u x w^r x 



(25) 



(26) 



(27) 



is an inverse operator to the above symplectic operator Mik iff 

1. f3 ik are rotation coefficients of a corresponding conjugate curvilinear coordinate net 
determined by the Bianchi-Darboux-Lame system (3) and by the "anti- flatness" condition 



diP ki + d k p ik + ^2p im P km = 0, i^k; 



{21 



2. affinors w l ,s determine N commuting hydrodynamic type systems (structural flows) 



r\ a — w^r x , 



(29) 



where wl-, = H^/Hi, and H( a y are solutions of the linear ODE systems 
diS(a)k = PikH( a )i, i ^ k, diH^y + ' s ^p im H^ m = 0; 

3. e alS is a constant non- degenerate symmetric matrix such that 
Proof: can be obtained by the direct verification 



(30) 



IC S M S] = 5}, 



M is K sj = Sl 



dxdt 



Corollary: The Lagrangian 

S = /^E^M^f-Mr) 
determines an integrable hydrodynamic type system (6), where (cf. (18)) 
Hi = H- /3 \/3, qp = y^^H(/3)mH m , diqp = ipiH^i. 



(31) 



(32) 



It means that the Hamiltonian density h (see (7)) is determined by the linear system 



dilPk = PkiA, ^ fc, d i H i + ^2Pim H n 



(33) 
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where dih = ipiHf, the momentum density (cf. (15)) 



(34) 



h 

2 

can be obtained by the replacement — *■ Hi — > Hi in the above system, where 
dji = ^>iHi. If r — > t, — > tpi, Hi — > Hi in the above system (cf. (14)), the above 
Lagrangian determines the hydrodynamic type system (5), where dih = ^-H*. 
Remark: Since the Lagrangian density 

L = L(r,r x ,r t ) 

determines an energy-momentum tensor, whose components are conservation laws of the 
energy and the momentum, respectively 

J.?L\ _ ( k dL\ f dL \ _ ( _ dL\ 

then the hydrodynamic type system (6) determined by the Lagrangian (31) possesses 
these conservation laws too, which are given by (see (34)) 

d T h = \d x (J2 Hi) , \d T (E H l) = 9 * (E H m H m - h) . 

The first example of integrable hydrodynamic type systems associated with the above 
Lagrangian representation was found in [43]. Indeed, under the point transformation 
(from physical variables to the Riemann invariants) 

dr± = du ± y/f"'(p)dp 
the Lagrangian (23) reduces to (31), where the Lame coefficients are given by 

Hi = ± I , 

and the ideal gas dynamics (19) reduces to the diagonal form 

— ± y?(r + -r )\r x , 

where the function ip(z) is given in the implicit form 



Taking N arbitrary solutions tpf^ of the adjoint linear problem (4), the hydrodynamic 
type system (5) can be re-written in the Hamiltonian form (see (27) and cf. (12)) 

dh 

at=e^(w a p ) x D-\w]) x — , (35) 

where conservation law densities and corresponding fluxes of structural flows (29) 

a% = «(a)) a (36) 

are determined by their derivatives 



dia? = i>^Hi, duo" = 



Mi- 



ll 



5 Mixed bi-Hamiltonian structure 

Let us consider an arbitrary semi-Hamiltonian hydrodynamic type system (2). It means, 
that characteristic velocities u l (r) satisfy the integrability condition (see [56]) 

c) k v i 'V 1 

d ~k — * = dk ~r — v iz £j^ k - 

V K — V 1 V 3 — V 1 

Following S. Tsarev (see [56]) the Lame coefficients H k are given by 

«9 fc lntf t = -^-, i^k. (37) 
v — v % 

Following G. Darboux (see [10]) rotation coefficients (3 ik are given by (26). Suppose 
the rotation coefficients (3 ik associated with a given hydrodynamic type system satisfy 
the flatness (9) and anti-flatness (28) conditions, simultaneously. It means that this 
hydrodynamic type system has two Hamiltonian structures 

rj = {r i ,Hi} 1 = {r i ,H 2 } 2j 

where components of the local (Dubrovin-Novikov) Poisson structure 

{r\r 1 } 1 = A i] 5{x-x') 

are given by (16) 

tii tii tiitl k 
and components of the local symplectic structure 

M 3k {r\r% = 5)5{x-x') 

are given by (25) 

Ma = —^-D x —^-, M ik \i^ k = HiH k ( ^ — — ^ 

X X \ X X 

Thus, the recursion operator of the second order is a product of local Hamiltonian and 
symplectic operators of the first orders 

R k = A ij M jk , 

and the above hydrodynamic type system has infinitely many local Hamiltonian structures 
of all odd orders 

i = ^—1 = AUjfr k A ks — = A ij M jk A ks M sn A nm ^^- = ... (38) 
and infinitely many local symplectic structures of all odd orders 

M l3 r\ = -^f M l3 Ai k M kn r r t = , M l3 ti k M kn A ns M sm r? = ^-f ... 
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Theorem 2: The nonlinear PDE system 

(39) 

@iPik + dkPki + ^y^PmiPmk = 0, ®i@H + ®h flits + ^^flimflkm = 0, * 7^ ^ 

zs an integrable system, which is a consequence of compatibility conditions following from 
the over- determined linear PDE system 

d t H k = (3 ik Hi, i ± k, XHi = di ( diHi + J2^nH n J + J>™ f <9 m # m + ]T/3 mn #„ 

^ k, AV>i = 9< I 0<^j + J^flni^n J + J2^im ( dm1p m + ^flnm^n J 



Proof: follows from the above construction. 

Indeed, suppose some particular solution of the first linear system (4) is given. 
Then corresponding solution ipf^ of the adjoint linear system (4) can be obtained from 
(33) 

The new particular solution H± is given by (8) 



(0) 
m ' 



Thus, an iteration of the above substitutions 
H\ k) = d^f +1) = d^ k+1) + * = 0,1,2,... (40) 

leads to semi-infinite series of solutions of linear problems (4). It means, that correspond- 
ing commuting flows and conservation law densities can be found too 



t k TT XI 



Hi 

Since the iterations (40) are invertible (see (18) and (32)), another semi-infinite series 
can be found in quadratures 

zj(-k-i) trip) (-fc) a (- fc ) 1(13)9^'^ <9/l (_fc) - u(-fc) 

Let structural flows (36) be written via flat coordinates of the first local Hamiltonian 
structure of the Dubrovin-Novikov type (see (12)) 

dh 1 



13 



where 

did? = i)f ) H l , diW* = i)f H {l) i, H {l)i = diip {7)i +^2p mi ip^ )m , d { h 7 = ^ {l)i Hi. 

Then an integrable hierarchy of hydro dynamic type systems can be written in the bi- 
Hamiltonian form 

= f* D 'i&r = ^W^K).^- ■ ( 41 ) 

Remark: More general compatible bi-Hamiltonian structure recently was investigated 
in [38]. However, the nonlocal Hamiltonian operator (27) can be effectively constructed 
for any given flat metric namely in the above (flat-anti-flat) case. Moreover, just in this 
mixed corresponding recursion operator is local. 

Indeed, since 

where 

7 7,M Ma) 

then an integrable hierarchy of hydrodynamic type systems possesses the local Hamilto- 
nian structure of the third order given by 

a% = r^nj-^PJ^'IK—j-^-. (42) 

Remark: In N component case, an inverse matrix t/^ 7 (a, a x .) is a very complicated 
expression with respect to the first order derivatives. By this reason, the Riemann invari- 
ants r h (see (38)) are most appropriate coordinates (for instance, in the Whitham theory 
for an arbitrary genus higher than 0), except the case when the Riemann invariants cannot 
be found explicitly, for instance for a dispersionless limit of integrable dispersive systems 
(see e.g. [48]). 

Remark: One of very powerful approaches in an integrability of the nonlinear PDE 
system (3) was presented in [8]. Two linear PDE systems (4) should be related by the 
linear ODE system (with respect to each fixed index i) 

Hi= [ad? + ... + c n (r)h/v 

If n — 1, this is exactly already well-known flat case (8). The case considered in [8] can 
be restricted on a compatible couple of the above transformations 

H? +1) = ic m d? + ... + c {n)l (r)}4 k \ 4 k) = [c {0)i d? + ... + c (fl)i (r)]tf« 

generalizing the case considered in this Section. 

Remark: The above transformation of the second order (see (39) and (40)) 

if< fc+1 > = (a? + ...)^ fc) 
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is a first example of more general transformations 



H\ k+1) = (d™ + ...)H\ k) 

never studied except n — 1 (see below, Section 7). 

Remark: Thus, the first "puzzle" hidden in [42] is a consequence of the existence of 
a mixed bi-Hamiltonian structure discussed in this section. It means, that the sub-class 
of so-called "separable" Hamiltonians for two component hydrodynamic type systems 
introduced in [42] has pure differential-geometric meaning via the language of conjugate 
curvilinear coordinate nets in N component case. 

Remark: The nonlinear PDE system (39) is well-known in classical differential 
geometry (see e.g. [5]). See also recent investigation in [54], where (39) is called 
0(2N)/0(N) x O(N)— system (see also the generalized wave equation in [53]). This 
nonlinear PDE system (39) is a natural generalization of the Egorov case (see below Sec- 
tion 7): if rotation coefficients (3 ik are symmetric, then (39) reduces to well-known system 
describing Egorov flat metrics. 



6 Multi-Hamiltonian structures 

A comparison nonlinear and linear PDE systems ((3), (4), (8) and (9)) describing or- 
thogonal curvilinear coordinate nets (3 ik with corresponding nonlinear and linear PDE 
systems ((28) and (33)) describing anti-orthogonal curvilinear coordinate nets $ ik leads 
to the simple identity 

hk = Pki- (43) 
Definition: Corresponding curvilinear coordinate nets are said to be mirrored. 
Since a description of orthogonal curvilinear coordinate nets is equivalent to a de- 
scription of anti-orthogonal curvilinear coordinate nets, then, in fact, all known theorems 
in the theory of orthogonal curvilinear coordinate nets can be simply reformulated for 
anti-orthogonal case. 

Interpretation: Any integrable (semi-Hamiltonian) hydrodynamic type system has 
an infinite set of conservation laws and commuting flows. 

Let us write N arbitrary conservation laws for N — 1 arbitrary nontrivial commuting 
flows in the potential form 



\Vn J 



hi h 12 
hi h 2 2 

h N i hjy2 



h 2 N 



h 



NN 



d 



( h \ 

t 2 

\t N J 



(44) 



It means, that N — 1 commuting flows 



H(i 



(i)< 



k = 2,3,...,N 



(45) 



possess N conservation laws 



d t khn = d t ih ik , 



l,2,..,iV, 
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where 

djKk = ip(i)jH(k)j- 

Let us construct another set of N conservation laws and N — 1 commuting flows 
determined by the transposed matrix 



(46) 





( Z1 ) 




hi 


h 2 i ■ 


■ h N i 




( Xi \ 


d 


Z2 




hu 


h 2 2 ■ 


■ h>N2 


d 


X2 




\z N ) 




hiN 


h>2N ■ 


■ h^N 




\X N ) 



It means, that N — 1 commuting flows 



k = 2,3, N 



(47) 



possess N conservation laws 

d x kh u = d x ih ki , i = l,2,...,N. 
Corresponding rotation coefficients are given by (see (43)) 

diH(i)k _ ^(i)fe 



0i 



ik 







k 't ■ 



because the new Lame coefficients = i^mi (cf. (45) and (47)). 

Such families of hydrodynamic type systems are said to be mirrored. The afore- 
mentioned relationship between these hydrodynamic type systems (44) and (46) is well 
known in the projective differential geometry as the "duality condition" (see [1] and other 
references therein). 

One of the most interesting questions of the classical differential geometry which has 
appeared at studying of semi-Hamiltonian systems of hydrodynamical type is the descrip- 
tion of the surfaces admitting not trivial deformations with preservation of principal di- 
rections and principal curvatures. Then the number of essential parameters on which such 
deformations depend, is actually equal to number of various local Hamiltonian structures 
of a corresponding hydrodynamical type system [20]. Such local Hamiltonian structures 
are determined by differential- geometrical Poisson brackets of the first order (see [14]). In 
the same paper a bi-Hamiltonian structure of the system of averaged one-phase solutions 
of Korteweg de Vries equation has been considered. Later multi-Hamiltonian structures 
of systems of hydrodynamical type were studied in [44] , [49] . 

Example [22]: The Lame coefficients (37) for the hydrodynamic type system 



r x , i = l,2,...,N 



are given by 
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where /ij(V) are arbitrary functions. The rotation coefficients (26) are given by 

/4V)# fe 1/2 







ik 1/2/ as ffO i 



where we use the notation 



H- = ] 



_ r m\l/2 



(49) 



The substitution these rotation coefficients in the flatness condition (9) yields the following 
functional-differential system (see [56]) 



^i^iP^k + T^^'iP^k + ^k^kPki + T^^'kPki + f'mPmiPmk ~ 0; 



(50) 



m^i,k 



where we use the notation (3% = d^/H?. 

Lemma [22]: The above functional- differential equation has a polynomial solution 
parameterized by N + 1 arbitrary constants 

Hi = C (rY + C^rY' 1 + C 2 ( r y- 2 + ... + CV- 2 (r*) 2 + CW-ir* + C N . 

Remark [44]: These Riemann invariants r k are nothing but elliptic coordinates (see 
[10]). 

Theorem 3 [22]: Semi-Hamiltonian hydrodynamic type system cannot possess more 
than N + 1 local Hamiltonian structures of the Dubrovin-Novikov type. 

Thus, the above hydrodynamic type system (48) has a maximum of local Hamiltonian 
structures of the Dubrovin-Novikov type. 

Remark: Bi-Hamiltonian structures of the Dubrovin-Novikov type are described in 
the set of publications (see [13], [21], [39], [47]). Semi-Hamiltonian hydrodynamic type 
systems connected with mtz/iz-orthogonal curvilinear coordinate nets, where the number 
of corresponding local Hamiltonian structures of the Dubrovin-Novikov type is less than 
N + 1 is not described yet. 

Since orthogonal and anti-orthogonal cases are dual to each other, then following 
theorem is valid. 

Theorem 4: Semi-Hamiltonian hydrodynamic type system cannot possess more than 
N + 1 anti-flat nonlocal Hamiltonian structures (see (27), (28)). 
Example: The hydrodynamic type system 



z — 1,2, 



^=(!> m - 2r< ) t 
has N + 1 local Lagrangian representations (see (31)) 



N 



(51) 



S a 



„fc\n— 1 



2^ n 



) x 



rn^k 



dxdt, 



n 



l,2,...,iV+ 1. 
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Indeed, the Lame coefficients are given by (cf. (49)) 

H°= Y[(r l -r m )-^ 2 . 

It means, that corresponding rotation coefficients /3 ik are mirrored (see (43)) to the rota- 
tion coefficients (3 ik associated with elliptic coordinates (see (48)). 

7 The Egorov flat hydro dynamic type systems 

Let us consider the semi-Hamiltonian hydrodynamic type system (2) possessing the local 
Hamiltonian structure of the Dubrovin-Novikov type (12) and the couple of conservation 
laws 

a t = b x , b t = c x . (52) 

Definition: The above Hamiltonian hydrodynamic type system is said to be Egorov 
flat. 

Theorem 5: The Egorov flat hydrodynamic type system has infinitely many local 
Hamiltonian structures of all odd orders. 

Proof: If the semi-Hamiltonian hydrodynamic type system (2) has the couple of 
conservation laws (52), then (see [49]) a is a potential of the Egorov metric, i.e. 

d t a = Hi dib = HiH h d iC = Hf 

and rotation coefficients (3 ik are symmetric (see [16]). It means, that both linear PDE 
systems (4) coincide to each other (thus, we can identify ip k and H k below). Also the 
flatness (9) and anti-flatness (28) conditions coincide and reduce to 

SPik = 0, (53) 

where 5 = Y^d/dr m is a shift operator. Thus, if the Egorov hydrodynamic type system 
has a local Hamiltonian structure of the Dubrovin-Novikov type, then another nonlocal 
Hamiltonian structure associated with a local Lagrangian representation exists automat- 
ically; and vice versa. 

The main statement of this section: If rotation coefficients f3 ik are symmetric 
and depend on differences of Riemann invariants (see (53)) only, then corresponding 
hydrodynamic type system (2) has an infinite set of local Hamiltonian structures of all 
odd orders. 

Remark: In general case, an infinite set of Hamiltonian structures can be constructed 
starting with a couple of given Hamiltonian structures (one of them must be invertible) 
due to Magri's theorem (see [36]). In the Egorov case, we need just one local Hamiltonian 
structure of the Dubrovin-Novikov type or one local Lagrangian representation (31). 

Indeed, the first transformation (8) (associated with the flatness condition (9)) reduces 
to (see [56], [57]) 

H\ k) = 5H\ k+l) - (54) 
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the second transformation (33) (associated with the anti-flatness condition (28)) reduces 
to (cf. (40)) 

H (k-i) = 5H (k)_ 

Thus, the transformation of the second order (40) factorizes in a product of the above 
transformations of the first orders. The well-known theorem (in fact written by L. Bianchi 
in [5]) is a consequence of such a factorization: 

Theorem 6: The nonlinear PDE system (Bianchi-Lame-Darboux-Egorov) system 
(3), (53) is integrable. The corresponding linear PDE system (which is equivalent to N 
linear ODE systems for each fixed index i) is given by 

diH k = p ik Hi, t^k, XHi = SHi. (55) 

Proof: can be obtained by a computation of the compatibility conditions di(djH k ) = 
dj(diH k ) for i ^ j ^ k and di(d k H k ) = d k (diH k ) for % ^ k. 

Let us consider the consistency of the linear problem (4) with the symmetry operator 
S = Y,s k (r k )d/dr k . Under the scaling dR h (r h ) = dr k /s k (r k ) this symmetry operator S is 
equivalent to the shift operator 5. 

Theorem 7: The first linear problem (4) is compatible with the symmetry operator 
S; i.e. the compatibility conditions di(djH k ) = dj(diH k ) for i ^ j ^ k and di(d k H k ) = 
d k (diH k ) for i ^ k are valid, where 

\H t = (S + Cl (r l ))H t , 

if and only if 

S(3 lk = {j s[)(5 lk . 

The second linear problem (4) is compatible with the symmetry operator S; i.e. the 
compatibility conditions di(djip k ) = dj(diip k ) for i ^ j ^ k and di(d k ip k ) = d k (diip k ) for 
i k are valid, where 

A^ = (S + c 4 (r*M, 

if and only if 

SP ik = (c k -d- s' k )(3 lk . 

Proof: can be obtained by a straightforward computation. 
In the symmetric case (ft ik = (3 kj ) the above formulas reduce to 

XHi = (S + s'^)/2)H u Sf3 lk + (s> + s' k )(3 ik /2 = 0. (56) 

Theorem 8 [44], [49]: The semi-Hamiltonian Egorov hydrodynamic type system 
(2) has a local Hamiltonian structure of the Dubrovin-Novikov type iff the symmetry 
restrictions (56) are fulfilled (and vice versa). 

Proof: Indeed, if the Egorov hydrodynamic type system has a local Hamiltonian 
structure of the Dubrovin-Novikov type, then the flatness condition (see (50)) 

Sidi(3 ik + -s^ + s k d k f3 ki + ^s' k [3 kt + s m(3 mi (3 mk = (57) 
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reduces to (56). 

Thus, if the functional-differential system (56) has M independent solutions (M < iV+ 
2), then a corresponding Egorov hydrodynamic type system has also M local Hamiltonian 
structures of the Dubrovin-Novikov type as well as M nonlocal Hamiltonian structures 
associated with local Lagrangian representations (31). 

Examples: The ideal gas dynamics (19) is the Egorov hydrodynamic type system. 
Indeed, the Egorov pair of conservation laws (52) is given by 

p t = {pu) x , (pu) t = {pu 2 + p 2 f"(p) - 2pf'(p) + 2f(p)) x . 

Just one symmetry operator 5 = d/du is connected with sole local Hamiltonian structure 
of the Dubrovin-Novikov type and with sole nonlocal Hamiltonian structure associated 
with a local Lagrangian representation (31). Thus, only one infinite series of local Hamil- 
tonian structures of all odd orders is connected with the ideal gas dynamics (19). 

Whitham equations (see [44]). The averaged N phase solution of the Sinh-Gordon 
equation is the Egorov hydrodynamic type system possessing sole local Hamiltonian struc- 
ture of the Dubrovin-Novikov type. The averaged N phase solution of the Korteweg de 
Vries equation is the Egorov hydrodynamic type systems possessing two local Hamilto- 
nian structures of the Dubrovin-Novikov type. The averaged N phase solution of the 
nonlinear Schrodinger equation is the Egorov hydrodynamic type system possessing three 
local Hamiltonian structures of the Dubrovin-Novikov type. 

Thus, the first Egorov hydrodynamic type system possesses sole infinite series of local 
Hamiltonian structures of all odd orders; the second Egorov hydrodynamic type system 
possesses two infinite series of local Hamiltonian structures of all odd orders; the third 
Egorov hydrodynamic type system possesses three infinite series of local Hamiltonian 
structures of all odd orders. 

Exceptional linear-degenerate hydrodynamic type system 

J> m r m - £*5>™) rl i = 1,2, 

Since the Lame coefficients Hi = 1/Sr m (see (37)), the potential of the Egorov metric 
a = — 1/Sr m and the rotation coefficients (3 ik = — 1/Er m , then the functional-differential 
system (57) has N independent solutions. Indeed, the substitution rotation coefficients 
in (56) leads to the following identity 



r m 



This system has the solution Sk{r k ) = r k and N — 1 parametric solution Sk = const, where 
Ssfc = 0. Thus, the above hydrodynamic type system possesses N infinite series of local 
Hamiltonian structures of all odd orders. N parametric local Lagrangian representation 
is given by 

So,r k + Sk r\ 

Let us restrict our consideration on N — 1 commuting Egorov flat hydrodynamic type 
systems written via flat coordinates or, whose first conservation law is given by (see [45]) 



S 



dxdt. 



dt-tai = d t ia. 
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where a 13 = g^a^ (see (12)). 

Theorem 9 [13]: These set of hydrodynamic type systems 



is determined by a solution of the WDVV equation 



d 3 F „ d 3 F d 3 F „ d 3 F 

^9™ a _„ a _ ra _„ = a _ na _ ra _„ 9™- 



da a daPda~> da^da^da^ da a da<^da^ da^da^da^' 
Corollary: Conservation law densities satisfy the linear PDE system 

_ 7r? d 2 h d 3 F _ 7 „ d 2 h d 3 F 



da^daH daPdatda* da^da^ dcfldatdaP ' 
Taking into account (this is a consequence from (58) when (3=1) 

d 3 F 

9a(3 ~ daJdaJdaT 

the above linear PDE system reduces (for (3 = 1) to 



Tern 



d 2 h d 3 F d 2 h 



a 11 = 

da l daT da^da^da^ da^da^ 

Since this linear PDE system is consistent with the symmetry (shift) operator 5 = d/da l , 
then this formula reduces to iterations (in quadratures) of all conservation law densities 
(see [13]) 

d 2 h ( 3 k+1) dhf d 3 F 

p — TflV P I, _ n 1 9 n — 1 9 N 

da»da< 9 dan daPdatdar' V ' ' '"' ' 

where first N conservation law densities hp are (for instance) flat coordinates a@. Thus, 
the generating function of conservation law densities /i(a, A) is determined by the linear 
PDE system 

A-**- = an 9h (59) 

whose compatibility conditions (see [13]) lead to the aforementioned WDVV equation. 

Theorem 10: The above linear PDE system is equivalent to (55). 

Proof: A differentiation (58) with respect to Riemann invariants (45) leads to the 
identity (see [13]) 

8 3 F ^ H 

in 

da a da' 3 da^ ~ ^ H [1)m ' 

Expressing derivatives of h with respect to flat coordinates via derivatives of h with respect 
to Riemann invariants (see (7), (17) and (11), where we must identify -0^^ = #( 7 )i), (59) 
reduces to (55). 
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In this case the bi-Hamiltonian structure (41) reduces to (see [38]) 

a __ ¥n dh^ _ _ K _ av _^ ( d 2 F \ . / d 2 F \ gfeg+ij 

atk ~ 9 Vx ~d^~ 9 9 9 {d&dtf J \dZd& J a ~ddT> 

because the structural flows (29) (see also (36)) are given by (58) and (cf. (10) with 
(30)) g alS = s 136 . A corresponding local Hamiltonian structure of the third order (42) 
reduces to (cf. (20)) 

af k = g p7 D x U a/3 D x U^D x -^, 

where 



Thus, the above local Hamiltonian structure of the third order is N component gen- 
eralization of (20). 

Remark: Ya. Nutku and P. Olver explained the origin of the second order recursion 
operator for the ideal gas dynamics (19) as a product of two Sheftel-Teshukov first order 
recursion operators in [42]. E.V. Ferapontov has proved (mentioned in [56], see also [25]) 
that the existence of the Sheftel-Teshukov recursion operator (see [52], [55]) is equivalent 
to the existence of the symmetry operator S (see (54)) in appropriate Riemann invariants. 
Indeed, the eigenfunction problem (59) can be written in the form 

dh / d 2 F \ dh 



da a 9 x \da a daP) 



da a daP J x da~i' 

Thus, the recursion operator (translating the gradient dh^ j da a to the gradient dh^ k+l ^ / da 1 
is given by 

the inverse recursion operator 

G% = U^g lP D x 

is connected with the Sheftel-Teshukov recursion operator (translating corresponding 
commuting flows) written via flat coordinates (cf. (21)) 

R% = g^DsG'D- 1 ^ = D^un (60) 

Thus, this observation (made by Ya. Nutku and P. Olver) is valid for iV component 
Egorov Hamiltonian hydrodynamic type systems (cf. (42)). 
Example: Let us consider the couple of commuting flows 



x ■ 



b t = d x {c + z bb ), b y = d x z ab , 

Q d x z aby Cy d x z aa , 
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where the function z(a,b) satisfies the associativity equation (see [13]) 



_ 2 

Zaaa ^abb ^aab^bbb- 

These Egorov hydro dynamic type systems have only one local Hamiltonian structure 

n dht dh dh x 

a* = D x-^~, bt = D x —, c t = D x — , 
oc do oa 

n dh 2 , n dh2 n dh 2 

a y dc > °v m , c y u x ^ 

of the Dubrovin-Novikov type in general case, where hi = be + z b , h 2 = c 2 /2 + z a . 
These commuting flows can be written in the potential form 




a 


b 


c \ 




t) 


b 


c + z bb 


Zab 


l-l 




c 


Zab 


Zaa / 







(61) 



where 

dtt = Q x dx + Q t dt + Qydy. 
Then a new local Hamiltonian operator of the third order (cf. (20)) 



D. 
D x 
D r 



where the recursion operator R is a purely differential operator of the first order (cf. (21)) 

R = D x {W x y l 
according to (60), and the matrix W is given by (cf. (61)) 




W 



8 Second Nutku— Olver's "puzzle" 

Let us consider a particular case of the ideal gas dynamics (19) determined by the special 
choice f(p) = p 7 /7(7 — 1)(7 — 2). The polytropic gas dynamics (see [41]) 

Pt = {P u )x, u t = I "tt 



2 7-1, 

has three local Hamiltonian structures of the Dubrovin-Novikov type 



p\ = Mi 1 A?\ ^ 
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where the Hamiltonians are 

,2 



,7 



fnr p 
2 7(7-1) 



pudx, 



H 2 



pdx, Hi 



udx. 



Corresponding components AV of the above local Hamiltonian operators are given by 



At 



D x 
D x 



.4, 



pD x + D x p (7 - 2) D x u + u D a 
D x u + ( 1 -2)uD x p'~ 2 D x + D x p'- 2 



( 



A- ? 



upD x + D x u p 



1^-2)^ + ^-1 



D, r . 



¥ 2 + ^ 



L p 7-i 



up^ 2 D x + D x up^ 2 



Then the first recursion operator 



Ri 



(7 - 1) u + (7 - 2H. J D T - 1 2p + ^D" 1 
2p 7-2 + ( 7 _ 2 ) pt-*p x D- x (7 - l)u + ^D- 1 



is a ratio of first two local Hamiltonian structures, i.e. A 2 = RiA±. However A 3 ^ R\A 2 . 
It means, that 



iMa = 2(7 - 1) 



/ 



V 



2wpD z + (up),, 



~ 1 ~12 U 7-1 \ x 



[\{l-l)u 2 + 2^}D x \ 
+ [|(7-2)u 2 + ^ 



/ 



-(7-2) 



< !\,. l p x + P X D X l u x 



D x l Px + u x D x l u x 



u x D^u x + p x D^U 



•.T-2 



7-2 



xx I 7—2 



x 

7-2 



is a nonlocal Hamiltonian operator of the Ferapontov type (see [19]). A comparison of 
the above nonlocal Hamiltonian operator R\A 2 and the third local Hamiltonian operator 
A3 leads to the purely nonlocal Hamiltonian operator (27) 



7 - 1 . 1 - o 

2- A 3 i?iA 2 

7-2 3 7-2 



u p 

7-2 



-1 



7-2 



u p 

7-2 



1 

1 



W x D- x W x a, 



where the matrix W is given by the particular choice f(p) = p° 77(7 — 1)(7 — 2) in (22). 
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Taking into account the identity (see (27), (41), (42)) 



B = AxK^Ax, 

where B is a local Hamiltonian operator of the third order (20) and K is aforementioned 
nonlocal Hamiltonian operator connected with the local Lagrangian representation (23), 
we obtain an interesting algebraic relationship 

B = A X (2^^3 - ^ AsAr 1 ^ M (62) 

for local Hamiltonian operators found by Ya. Nutku and P. Olver in [42]. 
The above relationship 

k = 2^A 3 - -^—A 2 A^A 2 

7 — 2 7 — 2 

has a simple differential-geometric interpretation (see [19]). Any semi-Hamiltonian hy- 
drodynamic type system has an infinite series of nonlocal Hamiltonian structures of the 
Ferapontov type. These Hamiltonian structures are associated with surfaces with a flat 
normal bundle. A reconstruction of corresponding surfaces with a flat normal bundle 
is a very complicated problem (see [7]) in general. Let ga = H? be a diagonal met- 
ric of the first Hamiltonian structure (i.e. A\ is written via Riemann invariants r k , see 
(16)), then g^' = Hf jr % is a metric of the second Hamiltonian structure (i.e. A 2 ), and 
g^ = Hf/(r 1 ) 2 is a metric of the third Hamiltonian structure (i.e. A 3 and A 2 A~^ 1 A 2 ). 
Thus, the third metric is associated with two different submanifolds. It means that cor- 
responding the two-dimensional Riemann space with a flat normal bundle is embedded 
into a four dimensional Riemann space in the second (nonlocal) case. 

Remark: Ya. Nutku and P. Olver have found one local Hamiltonian operator of the 
third order. However, the polytropic gas dynamics has three local Hamiltonian structures 
and three local Lagrangian representations (31) too, because this Egorov hydrodynamic 
type systems possesses three symmetry operators (the shift operator 5 = £<9/<9r m , the 
scaling operator S = Y,d/dr m , the projective operator P = E(r m ) 2 <9/<9r m , see (56) and 
[2]), where each of them is connected with one local Hamiltonian structure and with one 
local Lagrangian representation, simultaneously. Thus, the polytropic gas dynamics has 
three local Hamiltonian structures of the third order. Corresponding flat coordinates are 
found in [43]. Then two other similar algebraic relationships (see (62)) can be found. 



9 Generalizations 

Theory of orthogonal curvilinear coordinate nets associated with a flat metric can be 
extended on surfaces with a flat normal bundle [19]. In such a case two linear PDE 
systems (4) are connected by the nonlocal transformation of the first order (cf. (8)) 

Hi = 8^ + YjPw*l>m + d,P W = 
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where e a p is a constant symmetric non-degenerate matrix, H\ a ' are particular solutions 
of (4), the flatness condition replaces by the more general Gauss-Codazzi equation 



P TJ («) ZT(^) 



Thus, an obvious generalization of an anti-flatness condition and a corresponding 
nonlocal transformation of the first order between two linear PDE systems (4) can be 
given by 

^ = diHi + YsPmiHm + e«^°y , ^P (/3) = (63) 

where 

<9A; + d k (3 ik + ^ A„A m = £ Q ^l a Vf ^ 

Example: The hydrodynamic type system (51) has iV + l local symplectic structures 
related by the recursion operator (cf. [22]) 



R 



( r 1 



o \ 




/ 


' X 


^ X 


^ X 


\ 




1 




' .1- 


^ X 


^ X 






~ 2 












r N ) 




V 


X 


r N 

X 


r N 

X 


J 



D-\ 



(2) 



Indeed, it is easy to verify that M\ k 
are local. However, M^ +2 ^ 



(3) 



(N+l) 



R^ is no longer local. Then one can check that 
corresponding relationships between linear PDE systems (4) associated with higher sym- 
plectic structures M- k +1+n ^ are given by (63). 

Some nonlocal Hamiltonian structures of the Ferapontov type can be inverted. How- 
ever, corresponding symplectic structures are nonlocal too (see [37]). Thus, a list of all 
possible local Hamiltonian structures associated with semi-Hamiltonian hydrodynamic 
type systems is exhausted in this paper. 

Remark: The paper [l] particularly is devoted to a connection of Temple's sub- 
class with a linear-degenerate sub-class of hydrodynamic type systems. For instance, the 
linear-degenerate hydrodynamic type system (see [23]) 



is mirrored to Temple's hydrodynamic type system 
whose commuting flow is well-known electrophoresis hydrodynamic type system (see [4]) 



n 



r. . 



The first hydrodynamic type system has an infinite series of nonlocal Hamiltonian struc- 
tures parameterized by iV arbitrary functions of a single variable (see [23]); the second 
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hydro dynamic type system (as well as the third one) has another infinite series of nonlocal 
Hamiltonian structures parameterized by N arbitrary functions of a single variable (see 
[19]). Thus, the construction presented in the Section 5 can be extended on a nonlocal 
(Ferapontov) case [19]. 

Remark: The mixed case (cf. (39)) 

di(3 ik + d k /3 ki + PmiPmk = SapH^Hjp, di(3 ki + 8 k f3 ik + PimPkm = z^fH 

was considered in [53] for special simple choice e a p = b~ a p and e a p = 0. Such a case 
is called the generalized Sinh-Gordon equation (may be better to use another notation: 
the generalized Bonnet equation, see a corresponding comment, for instance, in [59]) 
associated with the symmetric spaces SO(2N, l)/SO(N) x SO(N, 1) (see details in [53], 
[54]). Thus, the above nonlinear system is associated with an isometric immersion of 
iV— dimensional submanifolds into a pseudo-Euclidean space. 



10 Non-hydrodynamic type systems 

Any semi-Hamiltonian hydrodynamic type system (5) has an infinite set of commuting 
hydrodynamic type systems (6) parameterized by N arbitrary functions of a single variable 
(see [56], [57]). However, higher commuting flows 

= Mj;, Mm), U z = f2(u,U x ,U xx ,U xxx ),... (64) 

can be constructed in some cases (see [26], [52], [55]). If some hydrodynamic type system 
has the local Lagrangian representation (31), then its higher commuting flows (64) are 
determined by similar Lagrangians 



S = J ^^2Hl(r)^-h(r,r x ,r XXl r 



dxdt. 



Corresponding higher order conservation law densities also were investigated in [56], [58]). 

Theorem 11 [56]: A semi-Hamiltonian hydrodynamic type system has the conserva- 
tion law density 

X 

if and only if the extra condition 

^H 2 k d k v k = 

is fulfilled. 

Taking into account (5), the above constraint reduces to 



^^{HkdkHk — H k d k H k ) = 
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Corollary: If rotation coefficients (3 ik are symmetric and depend on differences of 
Riemann invariants (see (53)) only (the Egorov flat case), then the above constraint 
reduces to 



k ■ 



where (see (33)) ip k = 5Hk and ^> k = 5Hk- 

Lemma: The set of commuting hydrodynamic type systems (58) possesses an infinite 
series of higher conservation laws and higher commuting flows. 

Proof: The Lame coefficients ifffcv of the Egorov flat hydrodynamic type systems 
(58) (written in Riemann invariants (45)) depend on differences of Riemann invariants. 
Thus, the above constraint is fulfilled. 

This infinite series can be constructed iteratively (see (38)). The first commuting flow 
of the third order is given by the above Hamiltonian Hi = J h(r, r x )dx 

1 5ri 

All higher commuting flows and higher conservation law densities can be obtained due to 
the second Hamiltonian structure (see (38)) 



M *W> ~TT = M i3 A^ k M kn r^ = Mj j A 3 M kn A ns M sm r™ , 



Thus, this result also generalizes the corresponding construction presented in [42]. 

Remark: All these higher commuting flows (see (64)) are integrable "dispersive" 
systems with a rational dependence on first derivatives. They are integrable, because 
they possess infinitely many Hamiltonian structures, conservation laws and commuting 
flows. Their solutions should be considered elsewhere. 



11 Conclusion and outlook 

In this paper we were able to find answers on a set of important questions. 

• Lagrangian formulation for nonlocal Hamiltonian structures of hydrody- 
namic type systems associated with the anti-flatness condition is pre- 
sented. 

• Existence of two Hamiltonian structures associated with flatness and anti- 
flatness conditions implies an infinite set of local Hamiltonian structures 
of odd higher orders. 

• Description of multi-Hamiltonian structures associated with anti-flatness 
conditions is equivalent to the description of local multi-Hamiltonian 
structures of the Dubrovin— Novikov type. 

• Existence of one local Hamiltonian structure for the Egorov hydrody- 
namic type system implies the existence of infinitely many local Hamil- 
tonian structures of odd higher degrees. 
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• Non-local Hamiltonian structures associated with the anti-flatness con- 
dition are extended on an arbitrary "co-dimension". 

The problem of description of local Hamiltonian structures of odd orders was formu- 
lated by S.P. Tsarev (see the end of Dr. of Science Thesis). Below we formulate a most 
general conjecture, which should be investigated elsewhere. 

The theory of local differential-geometric Poisson brackets 

{y(x)X(a/)} = .t' ? iu.u,.u,,,.. ...)5(x - x'), (65) 

where 

A«(u, u x , u xx , ...) = g*(u)D2 + ^'(uK^r 1 + (c k J H< x + cL«)D n x ~ 2 + - (66) 
starts from the first order (see [14]) 

{u\x)y{x')} = (g ij {u)D x + 6j"(u)«*) S(x - x') 
connected with hydrodynamic type systems 

ul = vl(u)u k x , k = l,2,...,N. 

Later these Poisson brackets were generalized on nonlocal case (see [19], [37]) 

KOzOX'Or')} = [A^(u,u x ,u xx ,...)+J2e al3 w i a D x 1 w^ 5(x - x'), (67) 
where functions u l (x,t) simultaneously satisfy commuting flows 

u t<* = w l a (u,u x ,u xx , ...). 

Moreover, in general case the existence of such nonlocal Hamiltonian structure means 
integrability of corresponding PDE system. 

Just the Poisson brackets of the first order (see [19]) 

{«<(*), ui(x')} = (g ij (u)D x + b^(u)u k x + ^ e^w\ a)k ulD x l w\ p)m u™) S(x - x') (68) 

are connected with hydrodynamic type systems. 

This paper particularly deals with the local Poisson brackets (65), where the operator 
A ij (u,u x ,u xx , ...) is given by (38) (cf. (66)). 

The main conjecture of this paper is following: a most general Poisson bracket 
associated with an integrable hydrodynamic type system is given by (67), where the local 
part (cf. (66)) A* J (u, u x , u xx , ...) is determined by similar quasi-rational functions as in 
(38); the nonlocal part is exactly as in (68). 

Remark: Aforementioned differential-geometric Poisson brackets are invariant under 
an arbitrary point transformation u 1 = u l (u), while the new Poisson brackets (38) are 
not invariant. Quasi-rational function means rational with respect to higher derivatives. 
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